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We consider a Dark Matter candidate particle of spin 3/2 with neutrino-like Standard Model 
strength interactions. In the Majorana case, the particle can account for all of the Dark Matter 
for a range of masses between 70-160 GeV, depending on the strength of the Higgs couplings. The 
elementary spin-dependent cross section on the nucleon is calculated to be 1.7 x 10 -2 pb, and does 
not depend on the mass or any additional parameters. The Dirac case is excluded by the very 
large coherent cross section. The amplitude at the nuclear level is of the purely isovector form. We 
make detailed predictions for differential and total rates of scattering on a variety of nuclear targets 
of interest to the current direct detection experiments. For heavy targets the annual modulation 
amplitude is predicted to be negative, which may help to determine the mass of the WIMP if and 
when data become available. 

PACS numbers: 12.60.-i 



I. INTRODUCTION 

The combined MAXIMA-1 pQ, BOOMERANG 0, DASI [3] and COBE/DMR Cosmic Microwave 
Background (CMB) observations [¥ imply that the Universe is fiat [3] and that most of the matter 
in the Universe is Dark [BJ, i.e. exotic. These results have been confirmed and improved by the 
recent WMAP data [7]. Comb ining the data of these quite precise measurements one finds: 

n b = 0.0456 ±0.0015, n cDM = 0.228 ± 0.013, fi A = 0.726 ± 0.015 . 

Since any "invisible" non exotic component cannot possibly exceed 40% of the above n CDU [3], 
exotic (non baryonic) matter is required and there is room for cold dark matter candidates or 
WIMPs (Weakly Interacting Massive Particles). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the 
observed rotational curves, see e.g. the review [3], it is essential to directly detect such matter. 
Until dark matter is actually detected, we shall not be able to exclude the possibility that the 
rotation curves result from a modification of the laws of nature as we currently view them. This 
makes it imperative that we invest a maximum effort in attempting to directly detect dark matter in 
the laboratory. Furthermore such a direct detection will also unravel the nature of the constituents 
of dark matter. The possibility of such detection, however, depends on the nature of the dark matter 
constituents and their interactions. 

Since the WIMP's are expected to be extremely non relativistic, with average kinetic energy 
(T) ss 50 keV(m WIM p/100 GeV), they are not likely to excite the nucleus, even if they are quite 
massive m WIMP > 100 GeV. Therefore they can be directly detected mainly via the recoiling of a 
nucleus (A,Z) in elastic scattering. Computing the event rate for such a process can requires the 
following ingredients: 

i) An effective Lagrangian at the elementary particle (quark) level obtained in the framework of 
the prevailing particle theory. The most popular scenario is in the context of supersymmetry. Here 
the dark matter candidate is the LSP (Lightest Supersymmetric Particle) [TUlU6| . In this case the 
effective Lagrangian is constructed as described, e.g., in Refs. [TUIU8) . At least till supersymmetry 
is discovered in the laboratory (LHC) other approaches are also possible, such as, e.g., technibaryon 
p~9| . 120 ) . mirror matter [2T| [22] . Kaluza-Klein models with universal extra dimensions (23J [24] 



ii) A well defined procedure for transforming the amplitude thus obtained, using the previous effective 
Lagrangian, from the quark to the nucleon level. To achieve this one needs a quark model for the 
nucleon, see e.g. [T5I |2"5H2"7] . This step is particularly important in supersymmetry or other models 
dominated by a scalar interaction (intermediate Higgs etc), since, then, the elementary amplitude 
becomes proportional to the quark mass and the content of the nucleon in quarks other than u and 
d becomes very important. 

iii) knowledge of the relevant nuclear matrix elements 28, 29 , obtained with as reliable as possible 
many body nuclear wave functions, 

iv) knowledge of the WIMP density in our vicinity and its velocity distribution. 

From steps i) and ii) one obtains the nucleon cross sections. These can also be extracted from 
experimental event rates, if and when such data become available. Typically, experiments obtain 
upper limits on the event rates, so one can get exclusion plots on the nucleon cross sections as a 
function of the WIMP mass. The extracted cross sections are subject to uncertainty of the nuclear 
and astrophysical inputs from steps iii)-iv). 

In the standard nuclear recoil experiments, first proposed more than 30 years ago |30j . one has to 
face the problem that the reaction of interest does not have a characteristic feature to distinguish 
it from the background. So for the expected low counting rates the background is a formidable 
problem. Some special features of the WIMP-nuclear interaction can be exploited to reduce the 
background problems. Such are: 

i) the modulation effect [3"T] . [3"2"H4"0] ii) backward-forward asymmetry expected in directional exper- 
iments, i.e. experiments in which the direction of the recoiling nucleus is also observed [4l] |42fl53] . 

In connection with nuclear structure aspects, in a series of calculations, e.g. in J35J [Ml [55J and 
references therein, it has been shown that for the coherent contribution, due to the scalar interaction, 
the inclusion of the nuclear form factor is important, especially in the case of relatively heavy targets. 
They also showed that the nuclear spin cross sections are characterized by a single, i.e. essentially 
isospin independent, structure function and two static spin values, one for the proton and one for 
the neutron, which depend on the target. 

As we have already mentioned, an essential ingredient in direct WIMP detection is the WIMP 
density in our vicinity and, especially, the WIMP velocity distribution. Some of the calculations 
have considered various forms of phenomenological non symmetric velocity distributions 39, 45 , 46 
and some of them even more exotic dark matter flows like the late infall of dark matter into the 
galaxy, i.e caustic rings [56TI6T)] . dark matter orbiting the Sun [35] and Sagittarius dark matter |61j . 
We will employ here the standard Maxwell-Boltzmann (M-B) distribution for the WIMPs of our 
galaxy and we will not be concerned with other distributions [62-65]. 

In addition to computing the time averaged rates, these works studied the modulation effect. 
They showed that in the standard recoil experiments the modulation amplitude in the total rate 
may change sign for large reduced mass, i.e. heavy WIMPS and large A. This may be exploited as an 
additional discriminant to deduce the mass of the WIMP from the data, if and when they become 
available. There exists a plethora of direct dark matter experiments with the task of detecting 
WIMP event rates for a variety of targets, such as those employed in XENON10 [66], XENON100 
[57], XMASS [55J, ZEPLIN [BJ, PANDA-X [70], LUX [71], CDMS [7J], CoGENT [73], EDELWEISS 
[73], DAMA [751 [75], KIMS [77] and PICASSO [TSJITpJ. 

In the present paper we will employ these well-understood ingredients and consider another at- 
tractive candidate, namely a spin 3/2 particle suggested in j80j |5T], which has some interesting 
properties and satisfies all the required constraints. Our detailed calculations of relic abundance in 
Section [Hi] show that the candidate can account for Dark Matter in the mass range between 70-160 
GeV, which is mildly dependent on two additional parameters. This particle couples to the nucleon 
via Z-exchange. Therefore it can lead to large spin independent and spin dependent cross sections 
depending on whether it is a Dirac or Majorana particle. In both cases the required cross-sections 
and differential and total rates are calculated with only one unknown parameter, namely the mass 
of the particle. The former leads to coherence of all the neutrons in the target, and can be already 
excluded by presently available data. If, however, it is a Majorana fermion, the time component 
of the neutral current contributes only at order /3 2 ss 10 -6 and the coherence is suppressed. In 
this case the direct detection process is dominated by the isovector, purely spin cross section. For 
the nucleon we calculate this cross-section in Section |IV1 with essentially no model uncertainty and 
obtain the value <7^ m = 1.7 x 10~ 2 pb. We present the differential and total direct detection rates 



for a variety of target nuclei, for both a light and a heavy WIMP in Section VI We conclude that 
the model, while being subject to the usual astro and nuclear physics uncertainties but in principle 
no adjustable parameters, can serve as a benchmark for experiments with nuclei which are sensitive 
to spin, including those with Xenon, and that one can be optimistic that the required sensitivities 
will be achieved in the near future. 



II. A SPIN 3/2 PARTICLE AS A DARK MATTER CANDIDATE. 

As we have mentioned in the introduction, the direct detection of dark matter is central to physics 
and cosmology. Among other things, it will establish the nature of the dark matter constituents 
and will discriminate among the various particle models. We will consider here a neutral spin 3/2 
particle [501 [HI] , which was proposed as part of an extension of the electroweak sector of the SM 
with gauge bosons of spin 2 and matter of spin 3/2. 

It is assumed that a minimum extension of the standard model (SM) takes place so that the left 
handed component is put in a isodoublet of the SU(2)i and the right handed partners in isosinglets, 

X L =( X ± > \ withF = -l, x« withr = 0, Z" withy = -2. (1) 

These spin 3/2 particles should be described by a Lorentz structure transforming as a vector-spinor; 
for brevity we have suppressed these indices. Thus, except for the different spin, it has the quantum 
numbers of a lepton doublet: x is neutral and I has unit negative charge according to the hypercharge 
assignments above. 

This particle possesses Yukawa couplings like every other fermion: 

£x = fx(x L J~) (f )x R , £ x = f x {xj-) (^ + o)c ( 2 ) 

and it can acquire a Dirac mass like every other fermion, i.e. spontaneously when the isodoublet 
acquires a vacuum expectation value 

< = k ^ (3) 

The neutral component may also acquire a Majorana mass like the neutrino, in a see saw like manner. 
This can be achieved directly in the presence of an isotriplet A of Higgs scalars [82 , 83J whose charge 
decomposition is <5 , 5~ , 5°. 



Then this leads to the coupling: 

£a, x = h A, X {xJl) ( I— § 
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This after the isotriplet acquires a vacuum expectation value becomes 

£a, x = h Aa (xJ~ L 
yielding the neutrino Majorana mass: 



o -va/V2 \ ( r R c 
o o \ - Xr 



m™=h A , x ^l (4) 

Alternatively the Majorana mass matrix can be obtained assuming that the isotriplet A possesses 
a cubic coupling with two standard Higgs doublets [84H87] (see Fig. [T]). This effective cubic 
coupling yields a Majorana mass term through the vacuum expectation value of the isodoublet. The 
effective Majorana neutrino mass is obtained from Eq. Q via the substitution 




Figure 1: The tree level contribution to the neutrino mass mediated by an isotriplet scalar. 



where is the vacuum expectation value of the standard Higgs doublet, wa is the mass of 6°. 

The cubic coupling may be thought of arising out of a quartic coupling A in the presence of an 
isosinglet acqiring a vacuum expectation value vr, in which case [aa = Au/j. 

We note here that, if the isotriplet A possesses a cubic coupling //a with two standard Higgs 
doublets, then this mechanism could lead to WIMP nucleus scattering via Higgs exchange, with 
one of the Higgs acquiring a vacuum expectation value. When it comes to the hadronic vertex this 
contribution is analogous to the LSP (light supersymmetric particle of supersymmctry). In any 
case, even if the situation is otherwise favorable, i.e. /ia xu« m^, this mechanism is going to be 
suppressed compared to the Z-exchange, due to the smallness of the mass of the dominant quarks 
present in the nucleus. We will not elaborate further on this mechanism. 

If this model is to account for Dark Matter, the Yukawa couplings should be chosen such that the 
neutral particle is lighter than the charged partner. In this case one may ask whether the particle is 
stable with respect to decay to SM particles. Fortunately, there is a natural explanation for stability 
of the lightest BSM particle. According to the quantum number assignment, there may exist a 
vertex for the decay of the neutral x m to a W + and an electron e~ (and in the Majorana case also 
a W~ and a positron e + ). Owing to the Lorentz structure, such vertex cannot be of the minimal 
type, and in fact in the model described in [81] is completely prohibited even if any possible loop 
effective contributions are included. Nevertheless, one may write down a dimension-5 operator such 
as 

Lint = ^fe^ X u W W+ - q v W+) + h.c. (6) 

The interaction is therefore suppressed by some high mass scale M*, which may be identified with 
the Planck scale. In this case the decay width is estimated to be V sa am x {jj 1 ^) 2 ~ 1CP 45 GeV 
which is sufficient to make the particle long lived on the cosmological scale. 



III. RELIC ABUNDANCE 



Some years ago, Lee and Weinberg [55J obtained a lower bound on the mass of a stable neutrino- 
like Dirac massive particle from the condition that it should not over-close the universe. With 
the couplings precisely those of the SM neutrino, and modern value for dark matter abundance, 
^DAth 2 = 0.1099 ± 0.0062 one obtains the estimate M v w 10 GeV with the cross section for annihi- 
lation at rest into a light fermion 

G 2 F Ml 

a v = 2 ^ ~ lpb. (7) 

The whole scenario became disfavored as such particle was not seen experimentally, and moreover 
LEP limited the number of neutrinos lighter than w 45 GeV to just three. Further, it was con- 
vincingly demonstrated by Kainulainen et al |89j that although the annihilation crosssection drops 
somewhat after the Z peak, nevertheless the relic abundance of such a particle does not reach the 
value necessary to overdose the universe. The dominant annihilation channel was found by them 
to be into W + ,W~ pairs for particles more massive than about 100 GeV. Moreover, the crosssec- 
tion keeps growing for high values of the mass, in part because of a large and growing coupling to 



the Higgs. This feature was criticised in [5U] on very general grounds, arguing that the thermally 
averaged cross section should drop like a/M x in all channels. Nevertheless, the calculations of Kain- 
ulaincn et al show that it comes close to being able to account for the Dark Matter just before this 
problematic channel opens up, f2 <s=s 0.1 for M x ps 80GeV. 

It is also interesting to ask whether the picture is different if this particle is assumed to be 
Majorana. In this case the annihilation crosssection to light fcrmions is suppressed at rest due to 
the vanishing of the vector current for a majorana particle |91l I92| ; we find that the annihilation 
into W + ', W~ pairs is also vanishing at rest. The resulting abundance is somewhat higher than in 
the Dirac case but does not reach the value necessary to account for all of the Dark Matter above 
the Z peak. We have calculated the abundance for a range of parameters in the case of our spin 
3/2 candidate, and found that the abundance may equal the presently well- measured Dark Matter 
abundance for a range of masses, M x pa 70 — 110 GeV. 




(a) majorana, spin 1/2 (b) majorana, spin 3/2 



Figure 2: The relic abundance as a function of the mass of the candidate particle for some representative 
value of the charged lepton mass and no doublet-triplet mixing. Horizontal axis: mass of the DM particle 
in GeV. Vertical: abundance, QnAih 2 . The model was implemented and abundances were calculated using 

Micromegas [93] . 

The relic abundance depends mildly on the mass of the charged partner l~ (eq. [TJ. It must 
be heavier than our neutral candidate, and highest abundances are achieved when it is 2 to 4 
times heavier than the neutral particle: if it is any lighter then co-annihilation effects reduce the 
abundance, and if it is much heavier there is less cancellations due to the W + , W~ t-channel diagram, 
also reducing the abundance. In addition the precise value depends also on the amount of Higgs 
doublet-triplet mixing which determines the relative contribution to the mass of eq.Q vs eq (|5]). 

With the Higgs mass of around 125 GeV, the Z,h channel opens up above M x = 110 GeV and in 
fact dominates the annihilation crosssection at high values of M x . Fortunately, as of the writing of 
this paper, we now know the value of the Higgs mass, but the couplings of this particle have not 
been yet verified experimentally. At LHC, the four lepton decay channel of the Higgs probes the 
ZZh coupling. We have used the Standard Model couplings, but in some extensions of the SM the 
ZZh coupling is modified significantly; this may extend somewhat the allowed range of masses to 
160 GeV (see Fig. [3). 
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Figure 3: The abundance in the case that the Higgs coupling to Z is substantially smaller than the 

tree-level SM value. Otherwise as in Fig. [2| 



Figure 4: The X3/2 interaction with quarks in the nucleus is mediated by Z-exchange. 



IV. THE DIRECT DETECTION FORMALISM 



The main advantages of our WIMP candidate are: 

• It can interact with the quarks via exchange of Z, a relatively light particle with well known 
interactions (Fig. [4]). 

• It can lead to both coherent and spin dependent modes. 

The rate for the coherent mode is proportional to a^ h N 2 , i.e. to the square of the neutron 
number of the target. 

The spin dependent cross section is purely isovector, i.e. the spin dependent rate is proportional 
to cr^ m (f2p — fin) 2 , i-e. on one particle model parameter, cr^ ln . Cl p , Q n are the proton and 
neutron components of the nuclear spin matrix elements which are known for most targets of 
experimental interest. The elementary nucleon cross sections, cr^ h and cr^" 1 are going to be 
calculated in this work. 



In the standard weak interaction the neutral current contact interaction takes the form: 

H = ^ftJ?° (8) 
J!^ c is the hadronic current at the nucleon level, which in the non relativistic limit is of the form: 



11 



J p C = \ ((-1 + 4sin 2 W ), g A v) « (0, 9 -j<t) , for the proton (9) 
i.e. the time component almost vanishes. Also 

j* c = ((--, - > for thc neutron ( 10 ) 

Note that the space component of the neutron is opposite to that of the proton. 

The current j M is associated with the WIMP, given in [M], except that here we also have chiral 
interactions, 

r = x u r^ p (A v + l5 A A ) x p (11) 

where T^ p can be replaced by 7^ r\ vp in this particular case. The calculated currents can be cast 
into a form similar to that of the nucleon with the use of the standard quantum mechanical spin 
3/2 angular momentum operators £3/2, 



31= ( Ay, 3L3/2AA) (12) 

here Ay and A^ are the vector and axial couplings of the particle. 
At this point we encounter the following possibilities: 

1. The coherent, spin independent (SI) rate. 

We will first consider the case that the spin 3/2 particle is a Dirac particle. In this case, it is 



natural to take the neutral current to be of the purely left handed form, with Ay — —Aa = 1. 
Assuming that the WIMP is a Dirac type particle and the vector coupling does not vanish, 
since Ay = 1, the nuclear cross section is dominated by its coherent component and is mostly 
due to the neutron. The coherent part of the neutron cross section can be cast in the form 

< 1 = G |^ ~ E \M {M limi ^M fl m f )\ 2 (13) 

Mi ,Mf ,mi,rrif 

where Mi,rrii are the oncoming spin projections of the oncoming 3/2 WIMP and nucleon 
respectively, while Mf,m/ are the corresponding quantities for the outgoing particles. M. is 
the invariant amplitude. The last factor of 1/8 arises from the averaging over the initial spin 
states of the WIMP and the nucleon. 

The spin independent amplitude is diagonal and takes the form: 

M (Mi, mi -> M f ,m 2 ) = —- S Mz ,M f 5 mi , mf (14) 



Averaging over the nucleon spins we get: 



2 



l E \M (Mi,mi^ M f ,m f ) | 2 = J 4x 2 (\) =1/4 (15) 



.oh G F m N 1 G F m N 

- -^r 4 = -^r (16) 

Thus in our model: 

cr^ h fa 1.66 x 10~ 3 pb (17) 

2. The spin dependent (SD) cross section. 

Next we will consider the case that this particle is of Majorana nature, i.e. the particle coincides 
with its antiparticle. In this case the time component of the WIMP current vanishes in the 
non relativistic limit, since effectively Ay = 0. The cross section is then of the purely isovector 
form, and is identical for the neutron and the proton due to their equal axial coupling: 

„ Bpi „ _ G F m% ^ 1 ^ |M , ; u ,„ ( ^ M/jW/) |2 (lg) 



I 7T 



where M.\ (Mi,m,i — » M/,m/) is given in the table |l] A straightforward summation over the 
nucleon spins results in: 



1 ( . 1 , 1 , 1 „ 4\ 5 

- 4x-+4x h4x-+2x- = — 

8^ 4 36 3 9 J 12 

The nucleon cross section is then 

spin _ G 2 F m 2 N 2 5 

-^r 9A u (19) 

Thus in our model, and with ~ 1.24, 

erj" ( (Dirac) « 4.25 x 1CT 3 pb, crj n ( (Majorana) « 1.70 x 10" 2 pb (20) 

i.e. in the Dirac case the spin-dependent cross section is not of interest to us, since it is sub 
dominant for all targets of interest: the spin independent mechanism dominates due to coherence. 
The Majorana case in interesting since in this case there is no competing coherent mode. Note that 
in our model the spin dependent cross section does not depend on the heavy-quark content of the 



Table I: The various matrix elements M. of the axial current interaction. Mj, Mf stand for the 
spin projections of the initial and final spin 3/2 WIMP, while ra^ mj are corresponding spin 
projections of the nucleon. The nucleon isospin projection was taken to be +1/2. 
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nucleon. We note that this large cross section in the case of the Majorana WIMP is close to the 
current limits extracted from the recent KIMS data [55] and the earlier ZEPLIN data 

2 x 10 -2 KIMS CsI(T£), 1.8 x 1CT 2 ZEPPLIN-III collaboration 

for the dominant cross section relevant to their target. 

It is now straightforward to find the cross sections at the nuclear level. In the first case, the 
cross-section is coherent and scales quadratically with the number of the neutrons: 

Kt^ = N 2 a^F 2 ( q ) (2i) 

where F(q) is the nuclear form factor and in the second case it is pure spin and is of the isovector 
form: 

< P u "L r = I (fip - n nf o#" Fu(g) (22) 

where-Fii(g) is the spin response function, normalized so that it is unity at zero momentum transfer, 
and Op, fi„ are the static values of the nuclear proton and neutron spin components for the target 
nucleus, normalized so that they become equal to 3 for the free nucleon. 



V. THE FORMALISM FOR THE WIMP-NUCLEUS DIFFERENTIAL EVENT RATE 

The formalism adopted in this work is well known (see e.g. the recent reviews [18j HZ])- So we 
will briefly discuss its essential elements here. Before calculating the direct detection event rate, we 
will first deal with the WIMP velocity distribution. To this end we will follow the steps: 

• one starts with such distribution in the Galactic frame. 

• one transforms to the local coordinate system: 

y — > y + v s + 6 (sin ax — cos a cos 7?/ ± cos a sin jv s ) , y — — (23) 
with 7 w 7r/6. u s a unit vector in the Sun's direction of motion, x a unit vector radially out 



of the galaxy in our position and y = v s x x. The last term in the first expression of Eq. (23) 
corresponds to the motion of the Earth around the Sun with 5 being the ratio of modulus of 
the Earth's velocity around the Sun divided by the Sun's velocity around the center of the 
Galaxy, i.e. vq rs 220km/s and S « 0.135. The above formula assumes that the motion of 
both the Sun around the Galaxy and of the Earth around the Sun are uniformly circular. The 



exact orbits are, of course, more complicated [351135], but such deviations are not expected to 



significantly modify our results. In Eq. (23) a is the phase of the Earth (a = around June 
Snd) 1 . 

One integrates the velocity distribution over the angles and the result is multiplied by the 
velocity v due to the WIMP flux. 

The result is integrated from a minimum value v m i n to the maximum allowed velocity v max - 
In general, the escape velocity v esc in our galaxy is estimated to be in the range 550km/s< 
i>esc < 650km/s. In our calculations we assumed for the M-B distribution v esc s» 620km/s, 
even though the value[5H] of 550 km/s, which results from an analysis of data from the RAVE 
survey |100j . would have been more appropriate. The obtained results are not sensitive to 
this value. v m i n is a suitable parametrization in terms of the recoil energy and the target 
parameters, namely: 



Am p E R 

v mm = ^ (24) 

where Am p is the mass of the nucleus, \i r is the reduced mass of the WIMP-nucleus system 
and En is the energy transfer to the nucleus. 

With the above procedure one obtains the quantity g(v m i n ). For the M-B distribution in the local 
frame it is defined as follows: 

g{v mm ,v E {a)) = 1 I"™* e ~^ + 2v.v E ( a) W E ( a ))/vZ v dv d ^ ^ = ^ (25) 
ve{o) is the velocity of the Earth, including the velocity of the Sun around the galaxy, Ve(o) — 



eo (v s + S (sin ax — cos a cos 73/ + cos a sin jv s )) (see Eq. (23)). The above upper cut off value in the 
M-B is usually put in by hand. Such a cut off comes in naturally, however, in the case of velocity 
distributions obtained from the halo WIMP mass density in the Eddington approach (62] , which, in 
certain models, resemble a M-B distribution [63j . Its precise value is not, however, important for 
the results of the present paper. 

Even though the differential rate is proportional [99 to g(v m i n , Ve(o)), for the benefit of the ex- 
perimentalists, we would like to make more explicit the dependence of the differential rate on each 
of the variables entering the expression 9(« mm ,VE(a)) and in particular to isolate the coefficient of 
cos a term, which will provide the interesting modulation amplitude and make the time dependence 
explicit. This approach will be even more useful, when one integrates the differential rate to obtain 
the total event rate. 

To this end, we will find it useful to expand g(v m i n) Ve(cx)) in powers of 6, keeping terms up to 
linear in <5 « 0.135. We found it convenient to express all velocities in units of the Sun's velocity vq 
to obtain: 

vog(v min ,v E (a)) = ^o(x) + ^i(x)cosa, x = (26) 

^q(x) represents the quantity relevant for the average rate and ^i(x), which is proportional to S, 
represents the effects of modulation. 

In the case of a M-B distribution these functions take the following form: 

^ (x) = i [erf(l -x) + eri{x + 1) + erfc(l - y OBC ) + erfc(y csc + 1) - 2] (27) 



One could, of course, make the time dependence of the rates due to the motion of the Earth more explicit by writing 
a « (6/5)tt (2(t/T) - 1), where t/T is the fraction of the year. 



where eri(x) and erfc(a;) are the error function and its complement, respectively, and y esc = v £SC /vq 
2.84. 

The differential event rate can be cast in the form: 
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cos a 



(29) 



where the first term represents the time averaged (non modulated) differential event rate, while 
the second gives the time dependent (modulated) one due to the motion of the Earth (see below). 
Furthermore for the coherent (here neutron coherent) in the case of a Dirac WIMP we get: 
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(30) 



with with /Lt r (fJLp) the WIMP-nucleus (nucleon) reduced mass and A is the nuclear mass number. 
m x is the WIMP mass, p(x) is the WIMP density in our vicinity, assumed to be 0.3 GeV cm -3 , and 
mt the mass of the target. Sometimes we will write the differential rate as: 
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V<v 2 > 



AT 2 ^"h 

(MA) N aN 



dt\ 
du ) 



r \ + H(ay/u) cos a) 



(31) 



In this formulation H(ay / u), the ratio of the modulated to the non modulated differential rate, gives 
the relative differential modulation amplitude, which is independent of the elementary cross section 
and the nuclear physics. 

Furthermore one can show that 



a 2 F 2 (u) * (a\/u), 



a 2 F 2 {u) Viiay/u) 



(32) 



The factor W2/3 is nothing but vq / U (v 2 ) since in Eq. ( 30 ) \/ (v 2 ) appears. In the above expressions 
a = (y/2/irbvo)^ 1 , vq the velocity of the sun around the center of the galaxy and b the nuclear 
harmonic oscillator size parameter characterizing the nuclear wave function, u is the energy transfer 
E R in dimcnsionlcss units given by 



E R 
QojA) 



Q (A) = [rripAb 2 ]- 1 = 40A" 4 / 3 MeV 



(33) 



and F(u) is the nuclear form factor. Note that the parameter a depends both on the WIMP mass, 
the target and the velocity distribution. Note also that for a given energy transfer E R the quantity 
u depends on A. The square of the form factor is exhibited in Fig. [5] 

For the axial current (spin induced) contribution, present for Dirac WIMPs and dominant for 
Majorana WIMPs one finds: 



dRo 



dE R 

dR 1 
dER 



Px m t 



m x Arrip 



Px m t 
m x Anip 



11 2 



Qo(A) 3 



dt 
du 

dh\ 
du J 



(34) 



(a) (b) 
E R -H>keV 

Figure 5: The square of the form factor entering the coherent scattering (a) and the spin response 
function in the isovector mode (b). With our normalization the other isospin modes are almost 
identical to the one shown. In both panels the solid, dotted and dashed curves correspond to the 
A=131, 73 and 19 systems respectively. The effect of the form factor for heavy targets becomes 
important even for an energy transfer of as low as 10 keV. 



dR 
with 



Px m t 
m x Am p 



< v z > 
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(l + H(ay/u) cos a) 



(35) 



: y|a 2 Fu(w) * (aV«), 



dh 
du 



y|a 2 ^n(u)*i(a^), (36) 



where Fix is the spin response function (the square of the spin form factor). The behavior of the 
spin responce function F\x for the isovector (isospin 1) channel is exhibited in Fig. [5j The other 
spin response functions Fbi and Fqq are, in our normalization, almost identical to the one shown. 
Note that, in the cases shown in Fig. [5] the spin response functions are not different from the square 
of the form factors entering the coherent mode. Thus, except for the scale of the event rates, the 
behavior of the coherent and the spin modes is almost identical. Anyway the quantity H(a^fu) 



appearing in Eqs (31 land (351 is the same 



Sometimes, the target has multiple components, either different atomic elements as in DAMA and 
KIMS or in the xenon experiments where several isotopes may contribute. In such cases the rates 
are additive, and the dimcnsionless recoil variable u must be normalized for each component: 



dR 
dE~n 



dR 



Er 



Qo{Ai) 



(37) 



where is the fraction of the component Ai in the target. 



Table II: The static spin matrix elements for various nuclei. For 3 He see Moulin, Mayet and 
Santos |103j . For the other light nuclei the calculations are from DIVARI [29]. For 73 Ge and 127 I 
the results presented are from Ressel et al [25J (*), the Finnish group et al jlOlj (**) and the 
Finnish team (+) [TU2] and Menendez (++) [TTJ4] . 



3 He 19 F 29 Si 23 Na 73 Ge 127 P 127 P* 129 Xe* 129 Xe+ 129 Xe++ 131 Xe* 131 Xe+ 131 Xe++ 

fio(O) 1.244 1.616 0.455 0.691 1.075 1.815 1.220 1.341 -0.609 1.174 -0.609 -0.325 -0.776 

fii(0) -1.527 1.675 -0.461 0.588 -1.003 1.105 1.230 -1.147 0.563 -1.105 0.567 0.320 0.679 

f2 p (0) -0.141 1.646 -0.003 0.640 0.036 1.460 1.225 0.097 -0.007 0.034 -0.023 -0.002 -0.023 

n n (0) 1-386 -0.030 0.459 0.051 1.040 0.355 -0.005 1.244 0.946 1.140 0.567 -0.323 -0.702 



Integrating the above differential rates we obtain the total rate including the time averaged rate 
and the relative modulation amplitude h for each mode given by: 



R mh = £* Jp- (tt) ^^N 2 a^t coh (1 + h coh co S a) , (38) 
m v Am p \HpJ 



i? api „ = ^ -P- (—) v^ 2 "^ I (n p - n n f off* t spin (1 + h spin cos a) . (39) 

m x Am p \HpJ o 

with 

t = I —du,h=- — du (40) 

JE th /Q (A) du 1 JE th /Q {A) du 

for each mode (spin and coherent). E t h{A) is the energy threshold imposed by the detector. 

The spin ME can be obtained in the context of a given nuclear model. Some such matrix elements 
of interest to the planned experiments are given in table [TlJ both in the isospin (fi ,f2i) and the 
proton neutron (f2 p ,f2 n ) bases. The shown results are obtained from DIVARI [3S], Ressel et al (*) 
[5S], the Finish group (**) [TOT] . (+) 



Using the expressions for nucleon cross sections, (38) and (39), we can obtain the total rates. 
These expressions contain the following parts: i) the parameters t and h, which contain the effect of 
the velocity distribution and the nuclear form factors ii) the elementary nucleon cross sections iii) 
the nuclear physics input (nuclear spin spin ME, see table [TT| . 

Let us briefly discuss the parameters t. These parameters depend on the reduced mass, velocity 
distribution and, for a heavy nucleus, on the nuclear form factors. In spite of this there is very 
little difference between the coherent and the spin contribution (see Fig. [6]) . The same thing applies 
in the case of the h factor. The target dependence will be inferred from the total rates computed 
below. 



VI. SOME RESULTS FOR MAJORANA WIMPS 



Before proceeding further we note that if the spin 3/2 particle were of the Dirac variety one 
would have a contribution arising from the time component of the neutral hadronic current. This 
essentially zero for the proton component, but it would have lead to a coherent contribution on the 
neutron. In this case the rate would be very large for heavy targets, proportional to iV 2 , which is 
excluded by the current data. We will, therefore, concentrate on the Majorana case, for which the 
coherent component will be proportional to /3 2 , (3 = v/c, and is therefore negligible in comparison 
to the spin-dependent component. 

We have seen that in our model the case of Dirac WIMPs is excluded. So we will concentrate on 
the Majorana spin 3/2 case. 
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Figure 6: The parameters t for the coherent (upper panels) and the spin (lower panels) in the case 
of the target A=131. On the right panels we exhibit these functions in the range of WIMP masses 
relevant to the present work, while on the left for comparison we show the low mass region. The 
curves correspond to threshold energies 0,1,2,4,7,10 keV, with the threshold energy increasing from 
top to bottom. We see that even for a heavy nucleus, for which the form factor becomes 
important, there is very little difference between the coherent and the spin t-factors. 



A. The differential event rates 

The differential event rates, perhaps the most interesting from an experimental point of view, 
depend on the WIMP mass. So we can only present them for some select masses. Considerations 
based on the relic abundance of the WIMP in this model lead to the conclusion that it has a mass 
between 80 and 200 GeV. This the WIMP mass range of interest to us. For illustration purposes, 
however, we have decided to present some results also for lighter WIMPs. Our results for the 
differential rates are exhibited in Figs[7}[9j. 

B. The total event rates 

We recall that the parameters h cob and h spin provide the ratio of the time varying amplitude divided 
by the time averaged rate. The relative modulated rate is given by h co h cos a and h sp i n cos a, with 
a the phase of the Earth. In the case of the spin we will present these functions in Figs [T0fT2"} 




5GeV 

lOCeV 

60GeV — -— - 
100 GeV - - - 
ISO GeV — — — 
250 GeV 



(a) 




E R ^keV 




E R ^keV 



Figure 7: On the top we show the time average differential rate ^Sr 1 



, and at the bottom panels 



we show the differential modulated amplitude as a function of the recoil energy Er in keV. 

These results correspond to the spin mode in the case of 131 Xe for a Majorana WIMP. On the left 
panels we exhibit curves for WIMP masses (5, 10, 60, 100, 180, 250) GeV, while on the right we 
exhibit WIMP masses relevant to our model, i.e. 60, 80, 180, 250 GeV. The relevant WIMP mass 

for each curve is explained on the top panels. 



As we have mentioned in the case of Majorana particles only the spin induced rates become 
relevant. Such results for the time averaged total rate are shown in Figs 13 - 15 . 



VII. DISCUSSION 



It is clear that, for the spin contribution, the target 19 F is the most favorable from an experimental 
point of view. This is due to the fact that the nuclear spin ME are large and the most reliably 
calculated. Such large rates favor the efforts of the PICASSO [75] ESI 1105) collaboration. As we 
have mentioned above our nucleon cross sections are consistent with the limits placed by KIMS data 
[95] and the earlier ZEPLIN data [96] . Thus we can conclude that relatively large rates predicted for 
these isotops are not ruled out so far. We have not yet seen an analysis of XENON 100 [57] in terms 
of spin induced cross sections. Our rates may seem a bit large, but in such an analysis it should be 
kept in mind that only 0.203 of the target is A=131. One has to also consider the A — 129 isotope. 
For the convenience of the experimental analysis we present the multiplicative factors needed to 
include to extend the data presented above in other isotopes and spin matrix elements in table |III| 
From table III we see that the rates may change by factors of 4, compared to the absolute rates 
presented in our figures. Even though the most recent calculation |104j obtains results which are in 
agreement with the earlier Ressel results, we should mention that these authors find a retardation of 
the isospin=l mode, which is the only one that contributes in our model, due to effects of long-range 
two-body currents in the nucleus [TU4] , [106] , This retardation, of about a factor of two in the rate, is 




E R -^keV 

Figure 8: The same as in Fig. Pmor the A=73 target. 
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>GeV 
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Figure 10: We show the relative modulation amplitude in the case of the spin for the target 131 Xe. 
Note the change in sign for heavy WIMPs, meaning that for negative sign the maximum occurs six 
months later than naively expected. As in Fig. |6j the curves correspond to threshold energies 
0,1,2,4,7,10 keV, ordered generally from top to bottom. 




essentially independent of the nuclear target. There is also the retardation of qa inside the nuclear 
medium (see e.g. ref. [T07J) , which may lead to another factor of 2. Taking this into consideration 
the rates predicted in the present work, which on the relatively large side, may be a factor of five or 
even 10 smaller due to uncertainties, which have nothing to do with particle model discussed here. 
In summary we have considered in this paper a viable and theoretically well motivated dark matter 
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Figure 12: The same as in Fig. 10 for the target F 





Table III: The multiplicative factor needed to go from the rates presented in this work to other 
possibilities of interest. These factors include the fraction of the isotope in question and the 

various spin matrix elements. 



Nuclear ME 


12a Xe(100%) 


ial Xe(100%) 


12a Xc(26.4%) , iai Xe(21.2%) 


Ressel 


1 


3.78 


1.21 


Finnish 


0.296 


2.61 


0.565 


Menendez 


1.32 


3.51 


1.67 



candidate, which satisfies all the experimental requirements. The results presented here may also 
apply to all models involving the WIMP interaction with quarks via the Z-exchange. 
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